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Abstract. Many materials quenched into their ordered phase undergo ageing and 
there show dynamical scaling. For any given dynamical exponent z, this can be 
extended to a new form of local scale-invariance which acts as a dynamical sym- 
metry. The scaling functions of the two-time correlation and response functions of 
ferromagnets with a non-conserved order parameter are determined. These results 
are in agreement with analytical and numerical studies of various models, especially 
the kinetic Glauber-Ising model in 2 and 3 dimensions. 



PACS: 05.70.Ln, 74.40.Gb, 64.60.Ht 

Ageing in its most general sense refers to the change of material properties 
as a function of time. In particular, physical ageing occurs when the under- 
lying microscopic processes are reversible while on the other hand, biological 
systems age because of irreversible chemical reactions going on within them. 
Historically, ageing phenomena were first observed in glassy systems, see pQ, 
but it is of interest to study them in systems without disorder. These should 
be conceptually simpler and therefore allow for a better understanding. In- 
sights gained this way may become useful for a later study of glassy systems. 

1 Phenomenology of ageing 

In describing the phenomenology of ageing system, we shall refer throughout 
to simple ferromagnets, see |2I3I4I5| for reviews. We consider systems which 
undergo a second-order equilibrium phase transition at a critical tempera- 
ture T c > and we shall assume throughout that the dynamics admits no 
macroscopic conservation law. Initially, the system is prepared in some initial 
state (typically one considers an initial temperature T ull = oo). The system 
is brought out of equilibrium by quenching it to a final temperature T <T C . 
Then T is fixed and the system's temporal evolution is studied. It turns out 
that the relaxation back to global equilibrium is very slow (e.g. algebraic in 
time) with a formally infinite relaxation time for all T < T c . 

Let 4>(t, r) denote the time- and space-dependent order parameter and 
consider the two-time correlation and (linear) response functions 



C(M;r) = Mt,r)0O»,O)) , R(t,s;r)= ^ff'^ 

on(s, 0) 



(1) 
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where h is the magnetic field conjugate to <f> and space-translation invari- 
ance was already implicitly assumed. The autocorrelation and autoreponse 
functions are given by C(t, s) = C(t, s; 0) and R(t, s) — R(t, s; 0) where t is 
referred to as observation time and s is called the waiting time. 
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Fig. 1. Two-time autocorrelator in the 2D Glauber-Ising model, with temperature 
T — 1.5 = 0.66T C on a 600 2 lattice and a disordered initial state. 



In figure^we show the autocorrelator C(t, s) of the 2D kinetic Ising model 
with Glauber dynamics after a quench to the final temperature T — 1.5. In 
the left panel, the dependence of C(t, s) on the time difference r = t — s is 
shown. Clearly, the autocorrelator depends on both t and s. For large values of 
s and r < s, the values of C(t, s) reach a quasistationary value C qs (r) ~ M c 2 q , 
where M eq is the equilibrium magnetization. In the regime r > s one observes 
an algebraic decay of C(t, s). Qualitatively similar behaviour is known from 
glassy systems and the simultaneous dependence of C(t, s) and/or R(t, s) on 
both t and s is the formal definition of ageing behaviour. The strong depen- 
dence of C(i, s) on the waiting time (which expresses the sensibility of the 
system's properties on its entire history) seems at first sight to lead to irrepro- 
ducible data and hence to prevent a theoretical understanding of the ageing 
phenomenon. Remarkably, Struik £Q observed in polymeric glasses submit- 
ted to mechanical stress that the linear responses of quite distinct materials 
could be mapped onto a single and universal master curve. We illustrate this 
here in the ferromagnetic Glauber-Ising model through the data collapse in 
the right panel of figure ^ Remarkably, a dynamical scaling holds although 
the equilibrium state need not be scale-invariant. 

On a more microscopic level, correlated domains of a linear size L(t) form. 
These are ordered if T < T c but do contain internal long-range fluctuations at 
criticality. In the first case, the system undergoes phase-ordering kinetics and 
in the second non- equilibrium critical dynamics. For sufficiently large times, 
the domain size scales with time as 



Lit) - t x ' z 



(2) 
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where z is the dynamical exponent. The slow relaxation to global equilibrium 
(although local equilibrium is rapidly achieved) comes about since for T < T c 
there are at least two distinct and competing equilibrium states. These states 
merge at T = T c . On each site r the local environment selects the local 
equilibrium state. 

As suggested from figure ^ one expects a scaling regime to occur when 

i *-§ > ^"micro ) S ^ > ^"micro ) t S ^micro (3) 

where r m i cro is some 'microscopic' time scale. We shall see later how important 
the third condition in is. If the conditions © hold, one expects [214] 

C(t,s) = S - b f c (t/s) , f c (y) ~ y- Xc/z ; y - oo (4a) 
R(t, s) = a- l - a f R (t/s) , f R (y) ~ y- x « /z ; y^ oo (4b) 

These scaling forms should hold for both T < T c and T — T c although 
the values of the exponents will in general be different in these two cases. 
Here Ac and \r are the autocorrelation [6J and autoresponse |7| exponents, 
respectively. They are independent of the equilibrium exponents and of z [Hj ■ 
It was taken for granted since a long time that Ac = A_r but examples to the 
contrary have recently been found for spatial long-range correlations in the 
initial data [7] and in the random-phase sine-Gordon model [§]. If Tj n i = oo, 
the inequality Ac = Xr > d/2 holds |TU]. 



Table 1. Values of the non-equilibrium exponents a, b and z for non- conserved fer- 
romagnets with T c > 0. The non-trivial critical-point value z c is model-dependent. 
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The values of the non-equilibrium exponents a and b apparently depend 
on properties of the equilibrium system as follows 1 1 1 j and are listed in 
table n together with those of z. We restrict to non-conserved ferromag- 
netic systems with T c > 0. If the equilibrium order parameter correlator 
C eq (r) ~ exp(— |r|/£) with a finite £, the system is said to be of class S and 
if C eq (r) ~ |t- | (rf— 2-i-tj) ^ ^ j g ^ ^ e Q £ c i ass £ At criticality, a system is 
always in class L, but if T < T c , systems such as the Glauber- Ising model are 
in class S, whereas the kinetic spherical model is in class L. For class S, the 
value of a comes from the well-accepted idea [1213] that the time-dependence 
of macroscopic averages comes from the motion of the domain walls. For class 
L, it follows from a hyperscaling argument 
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Having fixed the values of the critical exponents, we can state our main 
question: what can be said on the form of the universal scaling functions 
fc(y)i fn{y) ha a general, model-independent way ? 



2 Local scale- invariance 

Our starting point is the rich evidence, accumulated through many decades 
and reviewed in [2], in favour of dynamical scale-invariance in ageing phe- 
nomena. The order parameter field <f> = <p(t, r) scales 

(^(t,r) = b- x *(f>(b-'t,b- 1 r) (5) 

where b is a constant rescaling factor. We now ask whether eq. (JSJ can be 
sensibly generalized to general space-time dependent rescalings b = b(t,r) 
|14ll5llfi] . This ansatz can be motivated as follows. 

Example 1. Consider an equilibrium critical point in (1 + l)-dimensional 
space-time. Then z = 1 and let w = t+ir and w = t—ir. Any angle-perserving 
space-time transformation is conformal and is given by the analytic trans- 
formations w i — ^ f(w),w i— > f(w). A well-known result from field theory 
states that for short-ranged interactions, there is a Ward identity such 
that invariance under space- and time-translations, rotations and dilatations 
implies conformal invariance. Furthermore, basic quantities as the order pa- 
rameter are primary under the conformal group and transform as <fi(w, w) i— ► 
[f\w)f\w)) x *t 2 (j)(f(w)f(w)) (T^J. Hence n-point correlation functions and 
the values of the exponents can be found exactly from conformal sym- 
metry, see e.g. |17I19| for introductions. Here we merely need the projective 
conformal transformations f(w) — (aw + f3)/(^fw + 5) with aS — f3j = 1. 
Fields which transform covariantly under those are called quasiprimary |18|. 
The associated infinitesimal transformations are i n — —w n+1 d w , n £ {±1, 0}, 
and satisfy the Lie algebra [£„,^ n '] — (n — n')£ n+n >. 

Example 2. Let z = 2 and consider d space dimensions. The Schrodinger 
group Sch(d) is defined by [23j 

at + /3 Ur + vt + a 

1 1 ^ , r i ^ — — ; ad-Pj = l 6 

jt + o "ft + o 

where TZ € SO(d), a, v E R d and a, (3, 7, 5 € M. It is well-known that Sch(d) 
is the maximal kinematic group of the free Schrodinger equation Sip — with 
S = 2midt — d% [201 (that is, it maps any solution of Sip = to another solu- 
tion). There are many Schrodinger- invariant systems, e.g. non-relativistic free 
fields 21 or the Euler equations of fluid dynamics [221 • As in the conformal 
case, for local theories there is a Ward identity such that [23] 

space translation invariance 
scale invariance with z = 2 V Schrodinger invariance (7) 
Galilei invariance 
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We point out that Galilei invariance has to be required and for applications 
to ageing we note that time-translation invariance is not needed. Indeed, a 
non-trivial Galilei-invariance is only possible for a complex wave function ip. 
In applications to ageing, we shall identify below the 'complex conjugate' 
of the order parameter cf>* = cf> with the response field of non-equilibrium 
field theory 23 . We denote the Lie algebra of Sch(d) by scf) d . Specifically, 
sc §i = {^±1,0, ^±i/2; Mj} with the non-vanishing commutation relations 



where n,n' £ {±1,0} and m £ {±1/2}. 

Example 3. For a dynamical exponent z ^ 2, we construct infinitesimal 
generators of local scale transformations from the following requirements 16 
(for simplicity, set d = 1): (a) Transformations in time are t 1— > (at + (3) / (jt + 
S) with aS — f3j — I. (b) The generator for time-translations is JT_i = — dt 
and for dilatations Xq = —tdt — z~ 1 rd r —x/z, where x is the scaling dimension 
of the fields 0, <p on which the generators act. (c) Space-translation invariance 
is required, with generator — d r . Starting from these conditions, we can show 
by explicit construction that there exist generators X n , n £ {±1, 0}, and Y m , 
m = —1/z, 1 — l/z, . . . such that 



For generic values of z, it is sufficient to specify the 'special' generator \nj\ 



X X = -fdt-Ntrdr-Nxt-aSd^-pSd^-^-^-^r^lty 



from which all other generators can be recovered and where we wrote z = 2/N 
and 5, /3, 7 are free constants (further non-generic solutions exist for N = 1 
and TV = 2). For z = 2 we recover the Schrodinger Lie algebra sct) 1 . Now, the 
condition [X-y, Yjvya] — is only satisfied if either (I) j3 = 7 = which we call 
type I or else (II) a — which we call type II |16| . 

Definition: If a system is invariant under the generators of either type I or 
type II it is said to be locally scale-invariant of type I or type II, respectively. 

Local scale-invariance of type I can be used to describe strongly anisotropic 
equilibrium critical points. The application to Lifshitz points in 373 magnets 
with competing interactions is discussed in 24 16 . The generators of type II 
are suitable for applications to ageing phenomena and will be studied here. 
First, we note that the generators X n ,Y m form a kinematic symmetry of 
the linear differential equation Sip = where S = —z 2 ((3 + j)dt + |16| . 
Recently, systems of non-linear equations invariant under these generators 
with a = [3 — 7 = but extended to an infinite-dimensional symmetry 
1 1— > f(t) have been found [251 - Second, we consider the consequences for the 
scaling form of the response function R(t,s;r). To do this, we recall that 




(8) 




(9) 
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in the context of Martin-Siggia-Rose theory (see (2H|) a response function 
R(t, s) = (4>(t)cj)(s)) may be viewed as a correlator. If both (j) and 4> transform 
as quasiprimaries, the hypothesis of covariance of the autoresponse function 
leads to the two conditions X R — X\R = 0. Of course, ageing systems 
cannot be invariant under time-translations. From the explicit form of the 
generators given above these equations are easily solved and the result can 
be compared with the expected asymptotic behaviour (|4b(l . This leads to the 
general result |27llfi| 



R(t, s) = r o 0(t - s) ( - 



l+a-An/z 



(t- s y 



(ii) 



where tq is a normalization constant and the causality condition t > s is ex- 
plicitly included. Furthermore, the space-time response is given by R{t, s; r) = 
R(t, s)<P (r(t — s) -1 / 2 ) where <P(u) solves the equation 



% + z (/3 + 7) ud 2 - z + 2z{2 - z)7c)i- z ] ${u) 

In the special case z = 2, this reduces to ^3] 

M r 2 



= 



R(t, s; r) = R(t, s) exp 



2 t 



(12) 



(13) 



where M = j3 + 7 is constant. 

We point out that the derivation of the space-time response needs the 
assumption of Galilei- invariance (suitably generalized if z 7^ 2). In turn, the 
confirmation of the form l|13(l is a given system undergoing ageing provides 
evidence in favour of Galilei-invariance in that system. We shall next describe 
tests of (|11I13H in the Glauber-Ising model in d > 2 dimensions before we 
return to a fuller discussion of the physical origins of local scale-invariance. 



3 Numerical test in the Glauber-Ising model 

We wish to test the predictions (|11I13|) of local scale-invariance in the kinetic 
Glauber-Ising model, defined by the Hamiltonian H = a i a j where 

<7j = ±1. Based on a master equation, we use the heat-bath stochastic rule 

<7i(t + l) = ±l with probability i [l±tanh(/ii(i)/T)] (14) 

with the local field hi(i) — Xm(i) a n(t) and n(i) runs over the nearest neigh- 
bours of the site i. 

The response function is too noisy to be measured directly, therefore fol- 
lowing [2H| one may add a quenched spatially random magnetic field ±/i( ) be- 
tween the times t\ and £2 and measure the integrated response M(t, t\,t2) '■= 
/i(0) Jt^du R(t, u). Two schemes a widely used, namely the 'zero-field-cooling' 
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(ZFC) scheme, where t\ = s and t2 = t and the 'thermoremanent' (TRM) 
scheme, where t\ — and ti — s. However, in both schemes it is not possible 
to naively use the scaling form l|4bj) and integrate in order to obtain M. This 
comes about since in both cases some of the conditions (01 for the validity 
of this scaling form are violated. Taking this fact into account leads to the 
following results |11I13| : (a) the thermoremanent magnetization 

p(t,s):= [ duR(t,u) = r QS - a f M (t/ S )+r lS - x «/ z g M {t/ S ) (15) 



f M (y) = y-^'SFx ( 1 + a, ^ - a; ^ - a + 1; -) , g M (y) * ?T Ar/ * 
V z z yj 

where ro,i are normalization constants. The first term is as expected from 
naive scaling. In practice, a and Xr/z are often quite close and the size of 
the correction term may well be notable for T <T C (at T — T c , a and Xr/z 
are usually quite distinct); (b) the zero-field-cooled susceptibility 

X (t, s) J du R(t, u)= X o + s- A g{t/s) + O ( S - a ) (16) 

with a constant xo and some scaling function g. For systems of class S, we 
have A = a — k, where k measures the width w(t) ~ t K of the domain walls 
|13j . In the Glauber-Ising model, one has k = 1/4 in 2D and w(t) ~ Vint in 
3D 1211, while k — for d > 3. Consequently, the term of order s~ a coming 
from naive scaling is not even the dominant one in the long-time limit s — > oo 
and a simple phenomenological analysis of data of x(t, s) is likely to produce 
misleading results. For systems of class L, A = 0. 

Indeed, based on high-quality numerical MC data for x(t, s) in the 2D 
Glauber-Ising model and performing a straightforward scaling analysis ac- 
cording to xfa s) ~ s~ a but without taking the third condition @ for the va- 
lidity of scaling into account, it had been claimed that a = 1/4 in that model 

30 . However, that analysis is based on the identification A = a which cannot 
be maintained. Rather, for the 2D Glauber-Ising model, one has a — 1/2 and 
K = 1/4, reproducing A = 1/4 in agreement with the MC data. 



Table 2. Values of the autoresponse exponent Xr and of the parameters ro, ri and 
M in the Glauber-Ising model for an infinite-temperature initial state. 



d 


T 


Xr 


ro 


n 


M 


2 


1.5 


1.26 


1.76 ±0.03 


-1.84 ±0.03 


4.08 ± 0.04 


3 


3 


1.60 


0.10 ±0.01 


0.20 ±0.01 


4.22 ±0.05 



After these preparations, we can now present numerical Monte Carlo 
(MC) data and compare them with the predictions l|llll3fl . We consider the 



8 Malte Henkel 



thermoremanent magnetization p(t, s) and subtract off the leading finite-time 
correction according to l|15|) . For T < T c this leads to the parameter values 
collected in table see jHJ f° r details. Then the MC data both at T = T c 
and for T <T C for p(t, s) are in full agreement with (EJ 2ZEH, in both 2D 
and 3-D. Here we present a direct test of Galiliei-invariance by considering 
the space-time integrated response 



dp(t, s; p) 



= T du 



dr r^Rit, u; r) = r s d ' 2 - a p^ (t/s, p) (17) 



with an explicitly known expression for p^ following from l|13[) and the 
leading finite-time correction is already subtracted off Since all non- 
universal parameters were determined before and are listed in table [3 this 
comparison between simulation and local scale-invariance is parameter-free. 
The result in shown in figure [21 in 2D and we find a perfect agreement. A 
similar results holds in 3D |31j. 
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Fig. 2. Integrated space-time response of the 2D Glauber-Ising model. After 



This direct evidence in favour of Galilei-invariance in the phase-ordering 
kinetics of the Glauber-Ising model is all the more remarkable since the zero- 
temperature time-dependent Ginzburg-Landau equation (TDGL), which is 
usually thought to describe the same system (e.g. does not have this 
symmetry. Indeed a recent second-order result for R{t, s) does not agree with 
Ijlll) 32 (similar corrections also arise at T = T c [33]). On the other hand, 
Ac = 1 from the exact solution of the ID Glauber-Ising model at T — 0], 
while Ac — 0.6006 ... in the ID TDGL implying that these two models 
belong to distinct universality classes. 

Confirmations of (|1 1113(1 in exactly solvable models are reviewed in |16| . 
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4 Influence of noise 

We now wish to review the present state of theoretical arguments [28135] in 
order to understand from where the recent numerical evidence in favour of a 
larger dynamical symmetry than mere scale-invariance in ageing phenomena 
might come from. We shall do this here for phase-ordering kinetics. Then z = 
2 and we have to consider the Schrodinger group and Schrodinger-invariant 
systems. For simplicity, we often set d = 1. From the following discussion, 
the importance of Galilei- invariance will become clear, see also (0. 

A) Consider the free Schrodinger equation (2M.dt — d^)(f> = where 
A4 = \m is fixed. While an element of the Schrodinger group acts projectively 
(i.e. up to a known companion function 20 ) on the wave function (j), we can 
go over to a true representation by treating M. as an additional variable. 
Following |36| . we define a new coordinate £ and a new wave function ip by 

0(t,r) = -L [d(e- iM ^(C,t,r) (18) 

We denote time t as the zeroth coordinate and C as coordinate number —1. 

We inquire about the maximal kinematic group in this case |23| . Now, 
the projective phase factors can be absorbed into certain translations of the 
variable £ |23| . Furthermore, the free Schrodinger equation becomes 

(2id ( d t + d 2 r )^(C,t,r) =0 (19) 

In order to find the maximal kinematic symmetry of this equation, we recall 
that the three-dimensional Klein-Gordon equation J^„__ 1 d^d^^) = has 
the 3D conformal algebra conf 3 = so(4, 1) as maximal kinematic symmetry. 
By making the following change of variables 

C = (£o+i£-i)/2 , t=(-£o + i£-i)/2 , r = 6VV2 (20) 

and setting ip((, t, r) = !f r (^), the 31? Klein-Gordon equation reduces to (|19fl . 
Therefore, for variable masses M., the maximal kinematic symmetry alge- 
bra of the free Schrodinger equation in d dimensions is isomorphic to the 
conformal algebra conf d+2 and we have the inclusion of the complexified Lie 
algebras (sch d ) c C (conf d+2 ) c |23I37|. 

B) The Galilei-invariance of the free Schrodinger equation requires the 
existence of a formal 'complex conjugate' </>* of the order parameter </>. On the 
other hand, a common starting point in the description of ageing phenomena 
is a Langevin equation which may be turned into a field theory using the 
Martin-Siggia-Rose (MSR) formalism |26I17| and which involves besides 4> 
the response field (f>. If we identify <fi* = <fi and use l|18f) together with the 
assumption that ip is real to define the complex conjugate, then the causality 
condition that 



R(t,s;r) = (4>(t,r)4>(s,0)) - (0(i,r)^(a,O)> 



(21) 
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vanishes for t < s, follows naturally (and similarly for three-point response 
functions) [231 - Therefore, the calculation of response and of correlation func- 
tions from a dynamical symmetry should be done in the same way. 

C) So far, we have concentrated exclusively in applications of local scale- 
invariance to finding the form of response functions while the determination 
of correlation functions was not yet adressed. We shall do so now and consider 
the Langevin equation (with D -1 = 2M) |3S] 

d t0 = - D TT~ Dv{t)cp + ri (22) 
o<p 

where H. is the usual Ginzburg-Landau functional, v(t) is a time-dependent 
Lagrange multiplier which will be chosen to produce the constraint C(t, t) = 1 
and r\ is an uncorrelated gaussian noise describing the coupling to a heat 
bath such that (n) = and (r)(t, r)r)(t' ,r')) = 2DT5(t-t')S(r -r'). Another 
source of noise comes from the initial conditions and we shall always use an 
uncorrelated initial state such that 

C(0, 0; r) = (4,(0, r)<f>(0, 0)) = a S(r) (23) 

where clq is a constant. 

The MSR action of J23 reads S[<j>, <p] = S Q [<p, <p] + S b [<p, 4>] where 

S [<f>, 4>] = Jdtdr 4> + + Dv(t)^j (24a) 

S h [4>, 0] = -DT Jdtdr 4>{t, r) 2 - ^ J dr 0(0, r) 2 (24b) 

and we used (|23|l . see [32]. Here So describes the noiseless part of the action 
while the thermal and the initial noise are contained in S&. Finally, the po- 
tential v(t) can be absorbed into a gauge transformation; for example if 4>^ 
is a solution of the free Schrodinger equation, then <\> = (f>^k(t) solves the 
Schrodinger equation with the potential v(t) and where 

k(t) := exp (^-D J du v(u)j (25) 

The realization of the Schrodinger algebra with v(t) ^ is easily found |35j . 

We now assume in addition to dynamical scaling that TL is such that at 
temperature T = 0, the theory is Galilei-invariant |rS5| . This looks physically 
reasonable and we now explore some consequences of this hypothesis. We 
denote by (-)o an average carried out using only the noiseless part So of the 
action. The Bargman superselection rules state that (<fc- ■ ■ <j) 4* ' ' ' 0)o = if 

n m 

n =/= m. First, the response function is 

R(t, s; r) = (<j>(t, r)0(s, 0)) = (<f>(t, r)0(s, 0) exp (-S b [<p, 0]) ) 

= (#t, r)0(a,O)) =:-Ro(f, s;r) (26) 
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where in the last line the exponential was expanded and the Bargman super- 
selection rule was used. Here Rq is the noiseless response and the form (|1 11131) 
of Schrodinger-invariance is recovered if v(t) = (2M)(l + a — \n/2)t~ 1 |3"5] , 
In other words, under the stated hypothesis, the response function is inde- 
pendent of the noises. This is certainly in agreement with the explicit model 
calculations reviewed in section 3. 

Second, we now obtain the autocorrelation function. As before |35j 



where i?g (t, s,u;R) — (4>(t,y)(f>(s,y)4>(u, R + y) 2 )o- In contrast with the 
response function, the autocorrelation function contains only noisy terms and 
in fact vanishes in the absence of noise. By hypothesis, Schrodinger-invariance 

(3) 

holds for the noiseless theory and the three-point function ' is fixed up to 
a scaling function of a single variable [14135) . 

Working out the asymptotic behaviour of C(t, s) for y = t/s — > oo accord- 
ing to (|4a|) and comparing with the response function (|26|) . we find that for 
any coarsening system with a disordered initial state ]2ty and whose noiseless 
part is Schrddinger-invariant, the relation Ac = A_r holds true |35j . For the 
first time a general sufficient criterion for this exponent relation is found. 

The autocorrelator scaling function becomes for phase-ordering 



but the scaling function <P(w) is left undetermined by Schrodinger-invariance. 
If in addition we require that C(t,s) should be non-singular at t = s, the 
asymptotic behaviour <P(w) ~ u>~ Ac for w — > oo follows. Provided that form 
should hold true for all values of w, we would obtain approximately 



Indeed, this is found to be satisfied for several ageing spin systems with an 
underlying free-field theory [SSj. On the other hand, docs not hold true 
in the Glauber-Ising model. Work is presently being carried out in order to 
describe fc(y) m this model and will be reported elsewhere [35]. 
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